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Sunto. Si riconsidera il problema della rappresentazione di una mappa in una forma adatta 
all'applicazione dei metodi di forma normale. Si mostra che ricorrendo ai metodi delle serie di 
Lie e delle trasformate di Lie si pub costruire in modo diretto un algoritmo di normalizzazione. 
Si discute brevemente l'applicazione alia mappa di Schroder-Siegel e alia mappa standard di 
Chirikov, estendondole al caso di dimensione generica. 



Abstract. The problem of representing a class of maps in a form suited for application of 
normal form methods is revisited. It is shown that using the methods of Lie series and of Lie 
transform a normal form algorithm is constructed in a straightforward manner. The examples 
of the Scrhoder-Siegel map and of the Chirikov standard map are included, with extension to 
arbitrary dimension. 



Introduction 



Surface transformation as a tool for studying the flow of a system of differential equations 
have been introduced by Poincare ([20], Vol. Ill, ch. XXXIII) and deeply investigated 
by Birkhoff [2]. The corresponding method of Poincare sections has become classical, 
and has been widely used also for numerical explorations. A special interesting case is 
that of a periodic flow, for which the corresponding surface transformation is the time-T 
shift, T being the period. A natural question is whether a given map may be represented 
by the time-T flow of a differential system. Positive answers to this question exist for 
symplectic maps that are perturbations of integrable ones (see, e.g., [17], [18], [22], [13] 
and the references therein). 

In the present paper I will reconsider the problem of giving a suitable represen- 
tation of a class of maps with a method that is somehow connected to, but does not 
coincide with, the interpolation by a periodic flow. I will rather use the formalism of 
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Lie transforms, which allows us to extend to maps the techniques (e.g., calculation of 
normal forms) that are available for differential equations. 

The basic tools exploited here are not new: a short historical account is given at 
the beginning of sect. 2. However, most practical applications of Lie series and Lie 
transforms methods are related to differential equations. Typical subjects are numerical 
integration and coordinate transformations in perturbation theory. In the latter frame- 
work, in particular, these methods have proven to be very effective also in investigating 
the convergence or the asymptotic properties of perturbation series involving small di- 
visors (see, e.g., [7], [8], [9] and references therein) and in devising effective methods for 
perturbation expansions via algebraic manipulation on computers (see, e.g., [10] and 
references therein). 

It is a well known fact, however, that transporting the analytical methods of nor- 
mal form theory from differential equations to maps is not straightforward. A common 
remark is that the case of differential equations is easier to deal with, which justifies the 
attempts to interpolate a map with the Poincare section of a flow. 

A natural question is whether one can write a map in such a form that transporting 
the normal form methods that work fine for flows becomes straightforward. Answering 
this question in general is a major task, of course. However, if one considers a class 
of maps which are perturbation of integrable ones, then the question can be positively 
answered. 

One may consider as a basic example the Schroder-Siegel problem of iteration of 
analytic functions, extending it to many dimensions. However, the method developed 
here applies also to other cases, e.g., perturbations of integrable symplectic maps as 
considered in the papers quoted at the beginning of this section. As interesting models 
one may consider: the quadratic map investigated by Henon [14] and its generalization 
in higher dimension; the twist map of an invariant annulus investigated by Poincare [21], 
Birkhoff [1] and Moser [19]; the well know standard map. 

In all these cases the map may be represented as a composition of two maps: an 
integrable one (e.g., a linear one) and a near the identity perturbation. This is a trivial 
well known fact, of course. However, representing the integrable map as a Lie series and 
the perturbation as a Lie transform allows one to implement the normal form theory 
as a straightforward extension of the methods used for differential equations, provided 
a suitable formula for the composition of Lie transforms is available. This is what I'm 
going to illustrate. 

The paper is organized as follows. In sect. 2 a short account of the methods based on 
Lie series and Lie transforms is given, including the representation of a near the identity 
map and a composition formula. In sect. 3 the main proposition on the representation 
of a perturbation of an integrable map is proven. In sect. 4 it is shown how to construct 
a normal form algorithm for the map. The actual construction is worked out for two 
typical examples. A technical appendix follows. 

The theory is developed at a formal level. Some hints on quantitative applications 
to particular models are included at the end of sect. 4. 
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2. Basic tools 

The concept of Lie serie goes back to Sopus Lie. The use of Lie series in various problems 
has been widely investigated by Grobner in a series of papers after 1957 [11]. An accurate 
exposition with particular emphasis on applications to numerical integration can be 
found in Grobner's book [12]. The starting point, already found in Newton's work, is 
to express the solution of an holomorphic system of differential equations as a series 
expansion in time. The basic idea is to use the representation via power series as a one 
parameter near the identity map written in an explicit and useful form. 

Lie transform may be considered as a generalization of Lie series, in a sense that will 
be clarified later. Its usefulness as a tool in perturbation theory has been emphasized 
by Hori [16] and Deprit [3], who however paid attention in particular to Hamiltonian 
systems. The underlying idea of Deprit's work is to generate a one parameter family 
of near the identity coordinate transformation using the flow of a non autonomous 
system. It is quite common in the milieu of Celestial Mechanics to call Lie Transform the 
algorithms proposed by Hori and Deprit, reserving the name Lie series to the algorithm 
based on the flow of an autonomous system; thus I follow the tradition. 

Actually, several explicit algoritmhs for Lie transform have been proposed by many 
authors. A list of references may be found in Henrard's paper [15]. A similar algorithm 
for Hamiltonian systems has been introduced on a purely algebraic basis in [5]. In this 
paper I will follow the latter exposition, reformulating it for vector fields. A rigorous 
treatment may be found in [6], which has been prompted by the need of extending the 
contents of Grobner's book to a more general context. 

I recall here the definitions and the properties of Lie series and Lie transform 
working at a formal level, including what is needed in order to develop Lie methods 
for maps. I will omit most of the proofs, that can be found elsewhere. Furthermore, I 
will forget about the origin of Lie series as solutions of a system of differential equations, 
thus paying particular attention to the algebraic aspect. 

2.1 Lie derivatives 

Let V C C n be an open domain endowed with coordinates x = (x\, . . . ,x n ), and let 
X{x) = (Xl, . . . , X n ) be a complex holomorphic vector field on V. Let <f> x denote the 
time-t flow generated by X. 

The Lie derivative of a holomorphic function f(x) at the point x under the flow (f> x 
is the new function 



Similarly, the Lie derivative of a holomorphic vector field v is the new vector field 



(1) 




where ((f> x f)(x) = (/ o 4> x ) (x). 
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where {X, v} is the commutator between the vector fields X and v. 

The Lie derivative is a linear operator mapping the space of holomorphic functions 
(respectively holomorphic vector fields) into itself. It is also immediate to check also 
that the property 

L a X+/3Y = olL x + f3L Y 

holds true, where X, Y are vector fields and a, (3 are complex numbers. Further useful 
properties are the following. For two functions /, g the Leibniz rule applies, namely 

L x {fg) = fL x g + gL x f , L' x (fg) = £ (*) (L x f) {L 8 ~ j g) , s > 1 . 

3=0 

For two vector fields v, w one has 

L x {v,w} = {L x v,w} + {v,L x w} . 

Finally, denoting by [L x , Ly] = L x Ly — LyL x the commutator between the Lie deriva- 
tives with respect to the vector fieds X, Y one has 

[L X ,L Y ] = £{x,y} • 

The latter two properties are just different writings of Jacobi's identity for the com- 
mutator between vector fields, namely {X, {v, w}} + {v, {w, X}} + {w, {X, v}} = 0. In 
particular, if the vector fields X, Y do commute, namely if {X, Y } = then we have 

[L X ,Ly}=0. 

It will also be useful to write the explicit expression of the Lie derivatives in coor- 
dinates. For a function / one has 

n ft 

(4) 1> /• 

For a vector field one gets the expression of the commutator, namely 

where the l.h.s. is the j-th component of the vector field L x v. 
2.2 Lie series and Lie transform 

Let again V C C n be an open domain, and let X = (Xi, . . . , X n ) be a holomorphic 
vector field. The Lie series operator is defined as 

(5) exp(Lx) = J\ L x 

s>0 

A family of near the identity transformations depending on a parameter e may be 
constructed as 

(6) y = exp(eL x )x 
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or, in explicit form for the coordinates, 

e 2 

(7) y-j = exp(eL x )xj = xj + eXj(x) + —L x Xj{x) + ... , j = 1, . . . , n . 

The Lie transform is introduced as follows. Let X = {X 1 ,X 2 , . . .} be a sequence of 
holomorphic vector fields, that I will call the generating sequence; here the lower index 
labels the element of the sequence, not the component of the field in coordinates, that 
will be denoted, e.g., by Xij. The Lie transform operator is defined as 



(8) 



x 

s ' 



s>0 



where the sequence Ef of linear operators in recursively defined as 



(9) 



E, 



x 



i=i 



The superscript in E x is introduced in order to specify which sequence of vector fields 
is intended. However, I will remove it when unnecessary. By letting the sequence to have 
only one vector field different from zero, e.g., X = {0, . . . , 0, Xk, 0, . . .} it is easily seen 
that one gets T x = exp(Lx fc ). 

Writing e s X s in place of X s a one parameter family of near the identity transfor- 
mations may be defined as 



(10) 

i.e., in coordinates, 

(11) y 3 = Xj + eX^ix) + e 2 



y = T x x , 



-L Xl X hj (x) +X 2J (x) 



+ 



l,...,n 



Here I used the e expansion in order to make clear the connection with classical methods 
based on power expansion in a small parameter. However, in many cases it is convenient 
to just consider the vector field X s to be "small of order s" in some appropriate sense 
(e.g., using a norm) the order being implicit in the label of the field. Adding the powers 
of e the reader will easily check that every term in the definition of the operator Ef 
carries a factor e s , so that Ef is of order s. An equivalent precedure is to determine 
the order of E* as the sum of the indices of L x .Ef_-, which is s indeed. In the rest of 
the paper I will remove the parameter, unless it has a particular meaning. If the reader 
gets confused, he or she may just rewrite a formula by adding the powers of e, check 
that every ting is put in the correct order, and then set e = 1. 

The Lie series and Lie transform are linear operators acting on the space of holo- 
morphic functions and of holomorphic vector fields on the domain T>. They preserve 
products between functions and commutators between vector fields, i.e., if /, g are 
functions and v, w are vector fields then one has 



(12) 



Tx(fg) = T x f ■ T x g , T x {v,w} = {T x v,T x w} 
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Here, replacing Tx with exp(Lx) gives the corresponding property for Lie series. More- 
over both operators are invertible. The inverse of exp(Lx) is exp(— Lx), which is a 
natural fact if one recalls the origin of Lie series as a solution of an autonomous system 
of differential equations. The inverse of Tx takes the form 



yX 

s>0 

(13) 5 'J 

Gq = 1 , Gf = - -Gf_jLx n 



I come now to a remarkable property which justifies the usefulness of Lie methods 
in perturbation theory. I will adopt the name exchange theorem introduced by Grobner. 
Let / be a function and v be a vector field. Consider the near the identity transforma- 
tion (1 0) (or (6 ) for Lie series ) and denote by J the differential of the transformation (1 0) 
(or (6)), namely in coordinates, the jacobian matrix with elements Jj^ = ^f~. Then 
one has 

(14) f(y)\ =(T x f)(x) , rVy) =(Txv)(x). 

y y=Tx % 



This result should be interpreted as follows. On the l.h.s. of the equalities there is the 
transformed function (resp. vector field) calculated via the usual method of substitution 
of variables. The r.h.s. is the transformed function (resp. vector field) via the Lie trans- 
form, where the variables are renamed as x. The claim is that both operations give the 
same result. The remarkable fact is that if one uses expansions order by order, then the 
r.h.s. gives it in a straightforward way for the transformed function (resp. vector field) in 
terms of Lie derivatives, thus requiring only operations that are easily performed, e.g., 
via algebraic manipulation on computers. Obtaining the same result via substitution of 
variables is a definitely longer process, unless one stops the expansion at very low order. 
This claim may appear a little obscure, but it will turn out to be immediately evident 
if one writes the expansions as power series in a parameter. 

The proof of the identities (14) may be worked out using the algebraic properties of 
Tx - For, a holomorphic function (or vector field) may be expanded in power series of the 
variables, and by exploiting the linearity and the preservation of product the operator 
T x can be moved from the variables (the substitution) to the whole function. This 
justifies the name "exchange theorem" , since the symbol of the Lie transform operator 
is exchanged with the symbol of the function. 

2.3 Representation of a near the identity transformation 
Let now a near the identity transformation be given in the form 

(15) yj = Xj + eipij(x) + e 2 ip 2 ,j(x) + ... , j = 1, . . . , n , 

where <fiij(x), <p2,j(x), ■ ■ ■ are holomorphic functions. 
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Proposition 1: Any transformation of the form (15) may be written as a Lie transform 
of the coordinates y = Txx with the generating sequence 

^ k 

(16) X 1J = tp 1J , X rd = (p rJ ~22 -Lx k E r - k Xj , j = l,...,n, r > 2 . 

k=i 

In view of this proposition it may appear that Lie transform is more general and more 
attractive than Lie series, since the latter can not represent any near the identity trans- 
formation. However, a similar result may be obtained if one makes use of the composition 
of Lie series. 

Let again {Xi, X 2 , . . .} be a sequence of holomorphic vector fields. Consider the 
sequence of transformations {Sx\ Sx , S x \ ■ ■ • } recursively defined as 

(17) SP = 1, = exp(L Xr ) o . 
We may well consider in formal sense the limit 

(18) S x = ■ ■ ■ o exp(L Xr ) o . . . o exp(L Xa ) o exp(L Xl ) 

as an operator obtained by composition of Lie series. A noticeable fact, pointed out by 
Fasso [4], is that one has 

(19) S x = ■ E ? ■ 

s>0 

where : : denotes a reordering of in the following sense: a composition of r Lie 
derivatives Lx k Lx k2 • • • Lx kr is reordered as 

: L Xu Lxu o . . . o Lu : = Lx ,.,oLx o...oiv 

where cr(k) = {ai(k), . . . , a r (k)} is any permutation of k = {k\, . . . , k r } such that 
<7i (fc) > a 2 (k) > . . . > a r (k). E.g., : L Xl L X2 : = L Xa L Xl . 

Proposition 2: Any transformation of the form (15) may be represented via a compo- 
sition of Lie series of the form y = Sxx with a generating sequence X 2 , ■ . .} that 
can be explicitely determined with a recursive procedure. 

As the reader will notice, an explicit expression for the vector fields is missing in the 
statement. Such an expression may be produced exploiting (19), but it turns out to 
be quite useless, because it has a non recursive form. However, a recursive procedure 
for determining X\, X 2 , ■ ■ ■ may be easily constructed by trying the first steps. Setting 
Xij = ipij one has 

yj - exp(L Xl )xj = <P2,j - -L Xl X ld + ... 

where the dots denote terms at least of third order. Setting X 2 j = <p 2 j — \Lx x X\j one 
gets that yj — exp(Lx 2 ) o exp^Lx^Xj starts with terms at least of third order, which 
are used in order to determine X3, and so on. Such a procedure is easily implemented, 
e.g., via computer algebra, and it is in fact also the scheme of proof of the proposition. 
I add just a remark concerning the actual use of propositions 1 and 2. Apparently, it 
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seems unreasonable to do extra work in order to calculate the generating sequence of 
the transformation. However, if one wants to transform either a function or a vector 
field, pushing the calculation at high orders, then using the exchange theorem turns out 
to be definitely more effective than performing a substitution. 

2.4 Composition formulae 

I begin with a formula for the commutation between two Lie transforms. Let X, Y be 
generating sequences. Then one has 

(20) T x o T Y = T w o T x , 

with the generating sequence W = {TxY\, TxYi, ■ ■ •}• Here, the Lie trasform Tx may 
be replaced by exp(Lx) in case X is a vector field, and similarly for Y. 

Generally speaking the formula above does not seem very interesting because the 
vector field W turns out to be itself a series. However, it is useful in some case, typically 
when one is able to put W in a manageable form by explicitly calculating the sum of 
the series. It will be used later, in sect. 3.2. 

It is instead more interesting to to observe that since y = Tx ° Tyx is a near the 
identity transformation then in view of proposition 1 there exists a generating sequence 
Z such that y = T z x. The following proposition gives an explicit expression for Z. 

Proposition 3: Let X, Y be generating sequences. Then one has Tx ° Ty = Tz where 
Z is the generating sequence recursively defined as 

s-l . 

(21) Z 1 =X 1 + Y 1 , Z 8 = X s + Y s + J2 J -E?_ j Y j . 

3 = 1 S 

The proof requires a long sequence of tedious calculations that can hardly be found 
in previous papers. Thus I include it in appendix A. Composition formulas for any 
combination of Lie series and Lie transforms are easily obtained by suitably elaborating 
formula (21). E.g., if X, Y are vector fields then one has exp(Lx) ° exp(Ly) = Tw , 
where W = {W\, W2, . . .} is the generating sequence 

(22) W x = X + Y , W s = -L^Y . 

s . 

The latter formula reminds the well known Baker-Campbell-Hausdorff composition of 
exponentials. The difference is that the result is expressed as a Lie transform instead of 
an exponential. 

3. Representation of a map by a composition of Lie Transforms 

I shall consider two cases. The first one is a map in a neighbourhood of an equilibrium, 
which may be expanded in Taylor series as 

(23) z' = /\z + Vl (z) + v 2 (z) + . . . , zeC n 
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where A is a n x n complex matrix and v s (z) is a homogeneous polynomial af degree 
s + 1. The second example is a real analytic map 

<p' = <P + + e/ifo I) + e 2 f 2 (^ /) + ... 

/ / = / + £ Pl (^/) + £^ 2 (^/) + ... 

where (</?, i") G T n x C/, with C R m open, and £ is a small perturbation parameter. 
3.1 The unperturbed map 

By unperturbed map I mean here (as usual) either the linear part of (23) or the map (24) 
with e = 0, which is a Kronecker map on a family of invariant tori parameterized by 
the actions /. 

Let me start with the linear part z' = Az of the map (23). Let A = e A with a 
complex n x n matrix A. For the vector field Az one has L Az z = Az , . . . , L s Az z = A s z, 
and so also exp(L Az )^ = J2 s >o 7\ ^ z = eAz - Using the exchange theorem we may 
transform a function f(z') and a vector field v(z') as 

exp(L Az )f(z) =f(z') =/(Az), 

(25) z ' =e z 

exp(L Az )v(z) =e A v(z') =A 1 v(Az) . 

z'=exp(Z/Az)z 

The simplest case occurs when A = diag(Ai, . . . , X n ) is a diagonal matrix. 

I point out that here the exchange theorem is used only in order to represent a 
transformation of a function and of a vector field as the action of a Lie series operator. 
Of course, in this case everybody would perform the transformations via a direct sub- 
stitution as in the last member of the formulae above, since this is actually the sum of 
the Lie series in closed form and requires no further expansion. 

I come now to considering the map (24), that for e = writes 

(26) <p' = <p + , /' = I. 

Introducing the (n + m)-dimensional vector field O = (a; (J), 0) one immediately gets 

L a {<p,I)= (w(/),0) , L 8 v (u(I),0) = (0,0) for s > 1 . 

Thus, the map may be rewritten as (<//, I') = exp(Ln) (</?, I). Using again the exchange 
theorem as above we may transform a function f(ip, I) or a (n + m) -dimensional vector 
field v(<p, I) as 

exp(L n )/( V ,J) = /(<//,/') , (t =f(<p + u>(I),l) , 

e X p{L Q )v(<p,r)=r 1 v{<f/,I') =J-^(^ + u;(/),/) , 

(</?V )=exp(L n )(v,^) 

where J and its inverse J -1 are the jacobian block matrices 



In B ^ _ f\ n -B ^ g _ 



l m J ' V 'm / ' 1 dh 

\ s denoting the s x s identity matrix. 



1<?'<™, l<^<m 
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3.2 The full map 

The maps (23) or (24) considered in the last section are just interesting examples. 
In more general terms one can consider an unperturbed map x' = fo(x) which has 
some nice properties and can be represented as a Lie transform x' = exp(Lx) with an 
appropriately defined vector field X. In view of the particular role of the latter map I 
will denote R = exp(Lx). 

Then one considers a perturbed map 

(28) x' = Rx + f 1 (x) + f 2 (x) + ... 

where f s {x) is of order s in some reasonable sense. E.g., f s (x) should be replaced by 
v s (z) in (23), and by ef s ((p,I) in (24), with the obvious change of the symbols for the 
coordinates. 

Proposition 4: Consider the map (28) where R is a Lie series operator. Then there 
exist generating sequences of vector fields V(x) = |Vi(x), V 2 (x), . . .} and W(x) = 
[Wi(x),W 2 (x), . . .} with W s = RV S such that one has both 

(29) x = R o T v x and x' = T w oRx 

Proof. Using the linearity of the Lie series operator R rewrite the map (29) as 

x / = R(x + R- 1 (/ 1 + / 2 + ...)) . 

This is a representation of the map as the composition of two operations, namely 

<p(x) = x + R- 1 (A + f 2 + . . .) , x'= R0(x) . 

In view of proposition 1 one may determine the generating sequence V(x) = 
{Vi(x), V 2 {x)^ . . .} such that <p(x) = Tyx. Applying R to both members and using 
linearity one readily gets x' = R o Tyx , namely the first of (29). Using the iden- 
tity (20) of proposition 3, which clearly applies also to the Lie series operator R, one gets 
R o T v = T w o R with W as in the statement, which gives the second of (29). Q.E.D. 

3.3 Conjugating maps 

I come now to the following question. Let two maps 

(30) x' = T w oR x , y' = T z oRy 

be given, where R is an invertible Lie series operator and W = {Wi,W 2 , . . .}, Z = 
{Zi, Z 2 , . . .} are generating sequences. To find whether the maps are conjugated by a 
holomorphic near the identity transformation 

(31) y = x + <fi(x) + <f 2 (x) + . . . . 

Using the same operator R in both maps means only that the unperturbed maps are 
trivially conjugated. In view of propositions 1 and 2 it is natural to consider the trans- 
formation (31) as generated by either a Lie transform or a Lie series. 
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Proposition 5: Let X = {Xi , X 2 , . . .} be a generating sequence of the near the identity 
transformation y = Txx. Then the maps (30) are conjugated if 

(32) T w o T RX =T X oT z . 
More explicitly, the following relations must be satisfied: 

(33) DX 1 = Z 1 -W 1 , D = R - 1 

s-l . 

(34) DX S = Z s - W s + - {Ef-jZj - E^RXj) , s > 1 . 

3=1 S 

Proof. In the l.h.s. of the map y' = T z o Ry set y' = Txx' and then substitute 
x' = T\y oRx. This must be the same as substituting y = Txx in the r.h.s. of the map. 
Thus the identity 

T x x' 



= T z oRy 

x'—T w oR x 



y=T x x 



must be true, and by the exchange theorem this gives Tw ° R ° Tx = Tx ° Tz ° R. In 
view of (20) one has R o Tx = T R x ° R, so that (32) is readily found in view of the 
invertibility of R. Apply now proposition 3 to both members of the latter equality. By 
the first of (21) we get W 1 + RX 1 =X 1 + Z 1 , namely (33). For s > 1 we get 

s-l . s-l . 

•J TpW D ~V V I 'V I \ ^ 3 T7iX 



W s + RX S + J2 -Ef-j^Xj =X S + Z S + J2 -Ef-jZ 
s s 

3=1 3=1 



from which (34) readily follows. Q.E.D. 

A similar result for Lie series holds true if in (32) one replaces the Lie transform 
Tx with the infinite composition of Lie series Sx as given by (18). However, (34) must 
be restated in a more elaborated manner, proceeding step by step. Let me say that the 
maps (32) are conjugated up to order r in case there exists a finite generating sequence 
X = {Xi, . . . , X r } such that the transformation y = S^x makes the difference between 
the maps to be of order higher than r, i.e., 



-T z oRy 

x'=T w oRx 



y=S ( ^x 



0(r + l) , 



where is defined by (17). The maps are trivially conjugated up to order r if the 
generating sequences W, Z coincide up to order r, i.e., W\ = Z\, . . . , W r = Z r . For this 
implies T\yx — T z x = 0(r + 1). 

Proposition 6: Let the generating sequences of the maps (30) coincide up to order r — 1 
and let X r be a vector field of order r generating the near the identity transformation 
y = exp(Lx r x) . Then the maps are conjugated up to order r if 



(35) 



T w o exp(L R x r ) = exp(L Xr ) o T z 
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More explicitly, the following relations must be satisfied: 

(36) DX r = Z r -W r , D = R - 1 ; 

L(a-l)/rJ _ . 

(37) Z s - W s = -Ef_ r RX r - S —^V L x Z s-3r for s > r . 

The proof of (35) is a straightforward adaptation of that proposition 5 using a generating 
sequence X with all elements zero except X r . For, in this case one has Tx = exp(Lx r ), 
so that = if s is not a multiple of r. Removing from (34) all vanishing term, with 
some patience one obtains (37) . 

4. Normal form algorithm 

A standard and useful tool in perturbation theory is the construction of a normal form 
for either a map or a system of differential equations. This is a classical and widely 
investigated problem, so I will limit the discussion to indicating how the known methods 
may be revisited within the scheme of representation of maps presented in this paper. 
Constructing a normalization algorithm for maps is indeed an easy matter in view of 
the results of sect. 3.3. Actually, two different algorithms may be devised, the first one 
based on Lie transform, the second one based on composition of Lie series. 

4.1 A general formulation 

Let me start with the Lie transform. The key point is that formula (32) must be con- 
sidered as an equation for the generating sequence X, which must be so determined 
that the transformed sequence Z has some nice property that characterizes it as being 
in normal form. Look now at equations (34), which are just a rewriting of (32) order by 
order. The reader will immediately see that for every s > 1 one has to solve recursively 
an equation of the form 

(38) DX S + Z S = ^ S 

where \I/ S is known, since it is determined by Xi, . . . ,X 8 -i and Wi, . . . , W s -i, which 
are known. Thus, the problem is only that the prescription that Z has a normal form 
should be imposed so that eq. (38) may be solved for X s and Z s . This is the standard 
problem in normal form theory for both maps and differential equations. If the process 
can be worked out, at least formally, then the generating sequence X produces a coor- 
dinate transformation y = Txx, with inverse x = y, such that the map in the new 
coordinates writes y' = Tz o Ry. 

The algorithm based on composition of Lie series appears to be more elaborated, 
since it requires using proposition 6 as an iteration step. Precisely, one constructs an 
infinite sequence {W^} r >o of generating sequences, with = W, and a generating 
sequence X = {Xi, Xi-, ■ ■ .} such that for every r > the generating sequences 
and coincide up to order r — 1, and Z\ = W^~\ . . . , Z r = Wr are in normal 

form and do not change with the next iteration. Proposition 6 shows how to determine 
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X r and Wr = Z r by solving eq. (36) with W^ r ^ in place of W r . Then the whole 
generating sequence is constructed as given by (37), putting W^ r ~ 1S) in place of 

W and W^ r > in place of Z. The normal form is thus determined step by step. 

At first sight, the reader may think that this is a too complicated process. However, 
by implementing the algorithm using algebraic manipulation he or she will realize that 
there is no substantial increase of complexity with respect to the algorithm based on a 
single Lie transform, and that in some cases the composition of Lie series may even be 
more effective. 

4.2 Back to examples 

Let me illustrate how a normal form may be constructed for the examples (23) and (24). 
Actually, this means that I should explain how to characterize the normal form and how 
to solve eq. (38). 

Assume that the matrix A in (23) has a diagonal form, namely A = diag(Ai, . . . , A n ), 
and define Xj = e^ j+luJj . This is a n-dimensional version of the problem of iteration of 
analytic maps investigated by Schroder [23], who gave the formal solution for the case 
n = 1. Clearly one has 

R = exp(Lnx) , ft = diag(^i + iu u . . . , [i n + iu n ) , 

and applying the method of proposition 4 the generating sequence W may be determined 
so that W s is a homogeneous polynomial of degree s + 1. Following Schroder one tries 
to conjugate the map to its linear part. This means that the normal form of the map 
should be (' = AC, which in the Lie transform representation means that one wants 
Z = {0,0, . . .}, the null sequence. Thus, according to (38), the generating sequence is 
determined by solving for X s the equation 

(39) DX S = ^ S , D = exp(L nx ) -1 , 

where ^ s (xi, . . . , x n ) is a homogeneous polynomial of degree s+1. The relevant property 
is that D is diagonal on the basis of monomials x h ej = x\ x ■ . . .-x^ej, where (ei, . . . , e n ) 
is the canonical basis of C n . For, in view of the second of (25) one has 

Dx k e 3 = ( e <fc.A»+*">-A»i-*",- _ i) x k ej . 

Thus determining the vector field X s is an easy matter if none of the eigenvalues of D 
is zero. For, writing the homogeneous polynomial vector field as 

n 

3 = 1 |fc|=s + l 



the solution of (39) is readily found to be 



ti,h _ x k 



p{k,jj, J riuj) —jJij —iuj 

3 = 1 k 

Thus the map may be formally linearized if the nonresonance condition 
(40) e <fc,M+^)-M,-^ ^ i for k e Z ^ } \k\ > 1 and j = l,...,n 
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is fulfilled. In the case n = 1 considered by Schroder a resonance may occur only if A is 
a root of the unity. 

If the resonance condition is violated then a different definition of normal form 
must be introduced. Precisely, the space V s of homogeneous polynomials of any degree 
s splits into two complementary subspaces 

A/" 5 = D _1 {0} , n s = D(V S ) , 

namely the kernel of D and the image of V s through D. For, D maps V s onto itself, so 
that both M s and 1Z S are subspaces of V s , and D is diagonal. Then the operator D may 
be uniquely inverted on 1Z S . Thus eq. (38) may be solved by splitting \l/ s = tyg + 
with obvious meaning of the superscripts, and setting 

Z s = tff , X s = D~^f , X s e 1Z S . 

The latter condition makes the solution unique. The procedure thus described is a 
standard one in normal form theory. Different solutions may be considered, of course, 
depending on what one is looking for. 

The convergence of the tranformation to normal form in the non resonant case has 
been proved by Siegel [24] for the case n = 1 under the additional hypothesis that 
A satisfies a diophantine condition. The work of Siegel represents a milestone for the 
problem of convergencence of perturbation series with small divisors. 

Let me now come to the model (24). In view of the particular form of the map 
it is convenient to represent the vector fields separating, so to say, the ip component 

from the I component by writing f y 1 m P^ ace °f X, where X(<p, I) and Y(<p, I) are 

a n-dimensional and a m-dimensional vector function, respectively. The explicit form 
of the commutator is written as 



(41) 



x\ fv\\_( - vig) + S, {y,% - wijjf) 



F 'Iff 



Recall also that in this case whith the notation above one has 

D = exp(L n )-l, ^={ U( Q I) ) ■ 
In the r.h.s. of (38) we may expand * s (ip,I) in Fourier series as 

S foW^ J 

with known coefficients afc(7) and /3fc(J). Using a similar expansion for 

in view of (27) the action of the operator D is given by 

_ fE ke z4^ kMI)} - l)cfc(I)e*< fc '*> -E ke z^^ kMI)) Bd k (I)e^\ 
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Thus, one would be tempted to solve eq. (38) by setting 



(42) Z s 



a (I) 
A) (J) 



+ 



e i<fe,u>(/)>_l T ( e i(fe,^(/)>_ 1 )2 



Sfc€Z™\{0} 

V h _ . 

Z^A;6Z™\{0} e ^ fc ."( 7 )>-l 

However, some denominator could vanish at some point / G Q, the action's domain, 
or at least become very small. This is indeed the classical problem of small divisors in 
Celestial Mechanics, which was well known to, e.g., Lagrange and Laplace. 

Assume for a moment that no divisor actually vanishes. This is true, e.g., if u> e R n 
is a constant vector and the non resonance condition e l \ h ^) ^ 1 for k £ Z m \ {0} is 
satisfied. Then the generating sequence Z(I) is independent of the angle variables ip. 
Thus one has 

L z 3 (i)<P = Z S (I) , L r Zs ip = for r > 1 , 

and the map in normal form is written as 

if' = <p + u'(I) , l' = I + g'(I) 

where uj'(I) and g'(I) are determined via the normalization process. Thus the dynamics 
of the actions / is separated from that of the angles ip. 

A simpler form of the normalized map is found in case the map possesses some 



interesting symmetries. Let me give an example. Say that a vector field 

(+, — ) or of type (— , +), respectively, if it satisfies 

X(-<p,I)\ _ ( X(<p,I) \ _ fX(- V ,I)\ _ f-X(<p,I) 



is of type 



Y(-<p,I)J \-Y(^I)J \Y(-^I)J \Y(^I) 

With a little patience, looking at the explicit expression (41) of the commutator, one 
checks that the commutators obey the rules symbolically expressed by the table 



{•>•} 


(+,-) 


(-,+) 


(+,-) 


(-,+) 


(+,-) 


(-,+) 


(+,-) 


(-,+) 



Assume now that the map (24) satisfies the symmetry 

f 8 (-(p, I) = f 8 (<p, I) , g a (-<P, I) = -9s{v, I) ■ 

This means that the system is reversible. It is not difficult to check that in this case 
the generating sequence W(<p,I) is of type (+,—). Then, with a little more patience 
and using induction, one checks also that at every step of the normalization procedure 
one has that \l/ s is of type (+, — ), and so by solving eq. (38) one gets Z s of type (+, — ) 

and X s of type (—,+)• This implies that in (42) one has Z s = ^^q^^- -*- n turn this 

implies that the normalized map is written as 



ip + u'(I), I' = I 
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representing a Kronecker map on a family of invariant tori with angles oo'(I). 

All this is formal, of course. Making a rigorous statement, as is known, is a definitely 
more complicated matter. If the angles u(I) of the unperturbed map do depend on the 
actions I then the formal construction above is expected to fail due to the presence 
of zero divisors, and an analog of Poincare's theorem on non existence of holomorphic 
first integrals for Hamiltonian systems of differential equations applies. However one can 
prove, possibly with some extra condition on a>(I), that KAM theory applies, thus show- 
ing the existence of a big set on invariant tori carrying a Kronecker map with strongly 
non resonant angles. On the other hand, it should also be possible to prove a theorem 
of Nekhoroshev's type on exponential stability. If the angles u are constant, then some 
normal form may be constructed, possibly taking into account the resonances, but the 
series so constructed are expected to be divergent, unless one looks for a Kolmogorov's 
normal form on an invariant torus. However, all this matter goes behind the limits of 
the present note, which deals only with formal aspects. 



A. Proof of proposition 3 

Let me state a preliminary identity. If X is a generating sequence and V a vector field 
then one has 

s 

(43) EfL v = Y J L Efv Ef_ j . 

j=0 

The proof is worked out by induction, since the equality is trivial for s = and moreover 
for s = 1 it is just Jacobi's identity for commutators. Here is the complete calculation 
for s > 1. The notation is made simpler by writing Ej in place of Ef , since there is no 
confusion. 

s 

m 



J s — m—j 



E S L V — — L Xm E s - m Ly 

m=l 

s s — m 

= Yl 7 Lx - Yl L E 3 vE s -m-j 
m=l j=0 
s s—m 

= EE7 ( L Lx m E j V + L E . V L X S) E s 
m=l j=0 
s s—m 

= E E 1 j{ L Lx m E s _^ m vE J + L Ej vLx m E s - j - m ^ 

m=l j=0 

s-1 _ . s-j 

= Yl Yl J~( LL x m Es-^ m vEj + L Ej vLx m E a -j-r^j 

j=0 m=l J 

s-1 

Y ( L E B -jvEj + L EjV E s -j) 
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s . s — 1 ■ s 

= E -LejvEs-j + —^-L E . v E s _j = ^ L E . v E a _j . 

3=1 j=Q j=0 

Jacobi's identity for commutators written as LxL w — L w Lx = L^ xW is used in order 
to obtain the third equality. 

Coming to the proof of (21) of proposition 3, by definition of Lie transform one has 



(44) T X o T Y = £ E™ £ Ei P = E E Ei m X) E 



l>0 / K k>0 ' s>0m=0 

On the other hand for the generating sequence Z defined as in (21) one has 

(«) E ? = f +EjE jVvA • 

1=1 1=2 m=l 

Thus it is enough to check that 

s 

(46) E z s = E l E I-m , s>0. 

m=0 

I proceed by induction. For s = 0, 1 the equality is true. For s > 1 calculate 

s s — 1 



m=0 m=0 



V E x E Y = V S ~ m (E x E Y +E X E Y ) 



s — 1 s — m ^ 

E E { L x l E f- m -i E ™ +E X L Yl E Y _ m _^) 

m=0 1=1 



s — 1 s — m ^ / jn \ 

E E z[ L Xi E *-m-l E m + E L E?Y, E m-k E Y-m-l j 

m=0 1=1 ^ k=0 ' 
s s—l ^ 

= E E 7 { L Xl E f-m-l E m + L Yi E m E Y-m-l) 



S 

1=1 m=0 

s — 1 s — l m ^ 

+ E E E -LEXY^m-kK-m-l • 
2=1 m=l fc=l 

The identity (43) is used in order to obtain the third line. The first double sum in the 
latter expression is further elaborated as 

s—l s—l 



S £ / S—L S—L \ 

E ~ s ( E E f-l~m E m + L Yi E E m E I-l-m j 

1=1 ^ m=0 m=0 ' 



S I I 
= E ~ Lx i+ Y i E E ™ E I-l-m = E - L Xi+YtE?-l 
1 = 1 m=0 1=1 
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The latter expression coincides with the first sum in the r.h.s. of (45). The triple sum 
is further elaborated as 

s — l s — l ^ s — l — k s — l s—l ^ 

Y Y ~ s Le u y i Y E hEY-i- k -h = Y Y - L E^Yi^f-i- k 

1=1 k=l h=0 1=1 k=l 

s — l s ^ s m — l ^ 

= Y Y Z L E X ,Yi E lm = Y Y Z L E X ,Yi E s-m ■ 

f * f * § m — l 1 f * f * Q m — l b 

1=1 m=l+l m=2 1=1 

Here the induction hypothesis is used in the first step. The last expression coincides 
with the second double sum in the r.h.s. of (45). Thus, the right member of (46) equals 
the last member of (45), and this concludes the proof. 
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